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ABSTRACT 

Weak lensing surveys that can potentially place strong constraints on dark energy parameters can 
only do so if the source redshift means and error distributions are very well known. We investigate 
prospects for controlling errors in these quantities by exploiting their influence on the power spectra 
of the galaxies. Although, from the galaxy power spectra alone, sufficiently precise and simultaneous 
determination of redshift biases and variances is not possible, a strong consistency test is. Given 
the redshift error rms, galaxy power spectra can be used to determine the mean redshift of a group 
of galaxies to subpercent accuracy. Although galaxy power spectra cannot be used to determine 
the redshift error rms, they can be used to determine this rms divided by the Hubble parameter, a 
quantity that may be even more valuable for interpretation of cosmic shear data than the rms itself. 
We also show that galaxy power spectra, due to the baryonic acoustic oscillations, can potentially lead 
to constraints on dark energy that are competitive with those due to the cosmic shear power spectra 
from the same survey. 

Subject headings: cosmological parameters — large-scale structure of universe — method:statistical 



1. INTRODUCTION 



Prospects of cosm i c shear for cosmology Hul Il999t 
Hu fc Tegmarkl 1199ft IHu] 12001 IHutererl I2002t IHeavensI 



20031 IRefregieHTS IBenahed fc Van Waerhekel 
Tshak et al.1 120041 ISong fc Knral 120041 ITakada, fc .Tairl 
20041 iTakada fc Whitell2004l IKnox et al.ll2005D have in- 
spircd many on-going and future weak lensing surveys 
such as the Deep Lens Survey 1 , Canada France Hawaii 
Telescope Legacy Surv ey 2 , and the Large Synoptic Sur- 
vey Telescope 3 fLSST. IlVson et al.ll2003^ . 

The statistical properties of cosmic shear maps de- 
pend on the distances to the lenses and sources. With 
the rcdshifts of the source galaxies determined, the 
data can thus be used to constrain dark energy pa- 
rameters via their influence on the distance-redshift 
relati on l|Simpson fc Bridle! l2005t IKnox. Song, fc Tvsonl 
|2005|) . Because the source galaxies are faint and nu- 
merous their rcdshifts will be estimated from multi- 
band imaging data (photometric rcdshifts) rather than 
from spectroscopic data (spectroscopic redshifts). For 
redshift uncertainties not to degrade the constraints 
on dark energy parameters one must know the rms 
photomet ric redshift error and any biases very ac- 
curately llBernstein fc lair] 120041: iHuterer et all 120051: 



ti2005ft . 



Ushak fc Hiratall2005llMa. Hu. fc Hutere 

Photometric redshifts are traditionally calibrated with 
spectroscopy of a subset of the imaged galaxies. One 
needs a very large and fair sample though in order to 
achieve a sufficiently accurate calibration. For example, 
to reduce the uncertainties of the rms photometric red- 
shift error to 1%, a minimum of 5000 spectra are re- 
quired in each redshift bin and ideally for each spectral 
type of galaxy. This is achievable at low redshift, but 

Electronic address: zhan@physics.ucdavis.edu 
Electronic address: lknoxS!physics.ucdavis.edu 

1 Sec http://dls.physics.ucdavis.edu. 

2 See |http : //mm . cf ht .hawaii ■ edu/Scienc e/CFHLS/ . 

3 See |http : //mm. Isst ■ org. | 



it becomes prohibitively expensive toward high redshifts 
for the magnitude limits of the deepest planned surveys. 
Moreover, photometric redshift errors are not Gaussian; 
the distributions tend to have long tails. This means 
that the rms value may not be sufficient for characteriz- 
ing photometric errors and that one needs an even larger 
spectroscopic subsample to map the full error distribu- 
tion. Thus, it remains a challenge for very deep surveys 
to calibrate their photometric redshifts with sufficient ac- 
curacy. 

However this challenge is met, we expect that inde- 
pendent cross checks capable of revealing inconsisten- 
cies with their nominal error distributions, will be very 
valuable. The clustering properties of galaxies in red- 
shift space provide us with such an opportunity. Pho- 
tometric redshift errors, like peculiar velocities, alter 
the galaxy power spectrum in redshift space. Specifi- 
cally, they strongly suppress the power spectrum along 
the line of sight. The amount of suppression is expo- 
nent ially sensitive to the rm s photometric re dshift error 
(e.g.lSeo fc Eisensteinll2003l hereafter SE03; iZhan et alJ 
120051) . Thus" one can use the three-dimensional galaxy 
power spectrum to quantify photometric redshift errors 
and, consequently, prevent weak lensing statistics from 
severe degradation. 

The clustering properties of the galaxies are also 
dark energy probes in their own right. The fea- 
tures in the galaxy power spectrum (the broad- 
band shape and baryon acoustic oscilla tions) 

Peebles fc Yul Effil IBond fc Ef stath inil IT981 

[ Holtzmanl 1198ft lHu"fc Sugivama 119961 iCole et alJ 
120051 lEisenstein et all 120051 iHuetsil l2005|) can be used 
as CMB-calibratcd standard rulers for determining 
the angular-diamete r distance D(z) and con s train- 
ing dark energy jEisenstein. Hu. fc Tegmarkl fl_998t 
i Blake fc GlazebrooU 120031 IHu fe Haimanl 120031 iLinderl 
120031 SE03: IAngulo et a,l .II2005H . 

Several papers have studied the prospects for measur- 
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ing bary on acoustic oscil l ations f rom photometric surveys 
(SE03: IBla,ke fc Bridle) 120041 [HoTnev Jain, fc Takad^l 
l200llGlazebrook fe Blakel2005HLinderl2005|) . The main 
advantage of a photometric rcdshift survey, is the wide 
coverage, which reduces the sample variance error, and 
deep photometry, which leads to more galaxies and there- 
fore lower shot noise. 

The challenges for large photometric redshift surveys 
include redshift errors, dust extinction, g alaxy bias, 
redsh ift distortion and nonlinear evolution ijZhan et alJ 
IHH). These non-idealities do not produce strong os- 
cillating features in the power spectrum (similar to the 
baryon oscillations) and can be contro lled, despite some 
degradation to the measurement s l|Seo fc Eisensteinl 
I20M IWhiteil2T)?)Rt IZhan et al .112001 . 

Ideally, the clustering properties of the galaxies would 
be adequate to simultaneously constrain dark energy and 
the error distribution of the photometric rcdshift errors. 
As we show below though, such a self-calibration does not 
achieve a sufficiently accurate reconstruction of the pho- 
tometric redshift error distribution. However, the clus- 
tering properties of the galaxies can provide a valuable, 
though somewhat model-dependent, consistency test. 

In section 2 we present our method for forecasting con- 
straints on cosmological parameters as well as parame- 
ters of the redshift error distribution. In Section 3 we 
show how well the rms photometric rcdshift error can be 
estimated given the knowledge of the Hubble parameter 
H(z), and vice versa. We then demonstrate that mean 
rcdshifts can be accurately determined from the same 
data. In Section 4 we forecast constraints on dark energy 
from the broadband shape and baryon oscillations in the 
galaxies of the LSST survey. Since the constraints from 
the former alone are sub-dominant to those from the lat- 
ter (SE03), hereafter we only refer to baryon oscillations 
even though the broadband shape is always included in 
our analysis. In Section 5 we discuss these results and 
conclude. 

2. METHOD 

SE03 developed a two-stage Fisher matrix analysis to 
forecast errors on dark energy equation of state (EOS) 
parameters from baryon oscillations. We adopt their 
method with improvements and constrain photometric 
rcdshift errors at the same time. 

In the first stage of the SE03 analysis Hi and Di , where 
Hi = H(zi), Di = D(zi), and Zi is the mean redshift 
of the ith bin of galaxies, are treated as free parame- 
ters and their Fisher matrix is calculated. In the second 
stage, constraints on the matter density w m , the comov- 
ing angular diameter distance to the last scattering sur- 
face -Dcmb, Hi and Di are converted to constraints on 
dark energy parameters. The point of dividing the anal- 
ysis into two stages is to allow one to better understand 
how the dark energy constraints are arising. 

2.1. Observed Galaxy Power Spectrum 

Using a fiducial cosmological model as reference, we 
can write the observed galaxy power spectrum as 

x\W{c 2 a 2 k 2 /H 2 )\ 2 b 2 G 2 P(k) + P s , 



where is the linear redshift distortion parameter 
<IKaiserlll987fr . o~ z = o" 2 o(l + z) is the rms photometric 
redshift error, b is the galaxy clustering bias, G is the lin- 
ear growth function, P(k) is the matter power spectrum 
at z = 0, and P 8 is the shot noise. We have suppressed 
the argument z in functions D, D{, H, Hi, 0, cr z , b, and 
G for convenience. The true wavenumbers k± and fc|| are 
related to the fiducial ones by 

k± = k{±D { /D and k\\ = k { \\H/H { . (2) 

The window function, W, in Eq. is the Fourier 
transform in the radial direction of the real-space window 
function, W , which is the distribution of galaxies at true 
distance r given their estimated distance of f. With the 
assumption that this distribution has the form 

P(r\r)=W((r-f)/a r ) (3) 

the modulus of the Fourier transform only depends on 
fc||(T r = k\\ca z /H. If we further assumed a normal dis- 
tribution then \W\ 2 = e'^^V" 2 and Eq. (JTJ would be 
the same as in SE03. 

In all our calculations, for simplicity, we do assume a 
normal distribution for the redshift errors. But in gen- 
eral rcdshift errors arc significantly non-Gaussian. Typ- 
ically photometric redshift errors have long tails due to 
the fraction of redshift assignments that fail catastroph- 
ically. However, as can be seen from Eq. Q any non- 
Gaussianity has its impact on the observed galaxy power 
spectrum and, in principle, by assuming the isotropy 
of P(k), one can reconstruct the function W from the 
anisotropy of the observed galaxy power spectrum. Of 
course, with greater freedom allowed in the functional 
form of the redshift error distribution, the quality of the 
r econstruction will b e weakened. 

iScoccimarrol (|2004|) points out that the iKaiserl l|1987fi 
formula for the linear redshift distortion [used in Eq. 
may not be sufficiently accurate for the precision mea- 
surements we consider here. Fortunately, the effect 
of redshift distortion can be readily calibrated through 
./V-body simulations. Meanwhile, one can also quan- 
tify the window function W by assigning realistic rcd- 
shift errors to mock galaxies drawn from the simula- 
tions. Precisely calibrated matter power spectra are, 
in fact , requ i red for interpreting weak lensing data 
llWhitd 12001 IZhan fc Knoxl [200l iHuterer fc Takadal 
l2005HHagan. Ma. fc Kravtsovll2005|) . As such, calibrat- 
ing the effects of redshift distortion and photometric red- 
shift errors can be achieved from the same set of simula- 
tions without much difficulty. 

Although we treat them as constant across each red- 
shift bin, all the functions (D, Df, H, Hf, 0, a z , b, and 
G) in equation as well as the survey selection func- 
tion, can change considerably within a single redshift bin. 
As such, the observed galaxy power spectrum in a red- 
shift bin becomes a convolution of equation with an- 
other window function th at accounts for th e radial evo- 
lution (e.g., as treated in IZhan et alJ l2005|) . We avoid 
this complication for simplicity, although the appropriate 
treatment is straightforward as long as the true evolution 
within each bin is known, or can be parameterized. 

The sub-bin evolution we are most concerned with is 
that of the clustering bias. The variation will be caused 
in part by any redshift dependence of the selection func- 
tion. Even if one were selecting the same distribution of 



Photometric Redshift Consistency 



3 



galaxy types at each redshift, bias would change due to 
the intrinsic evolution of the population. Most impor- 
tantly, redshift errors might actually be correlated with 
bias fluctua tions since both can be caused by spectral 
type errors i|Padmanabhanll2005l) . We do not yet know 
the magnitude of these effects and the extent to which 
they can be controlled. 

2.2. Fisher Matrix 

For the stage I analysis, our complete parameteriza- 
tion of the observables includes more than just Hi and 
Di. There are four other parameters in each redshift bin: 
ln(G;6;), ln/3j, P s a and \a.a z ,i 1 where i refers to the ith 
redshift bin, and seven redshift-independent parameters. 
We do not include the redshift bias parameters, 6zi, at 
this stage because they have no influence on the observed 
power spectra except through Hi and Di, which are in- 
dependent parameters in stage I. They will be included 
in stage II. For brevity, we suppress the subscript i in 
what follows. 

The seven redshift-independent parameters are the 
matter density w m , the baryon density cjb, the redshift of 
rcionization z re i, the primordial helium mass fraction j/ p , 
the spectral tilt n s , the normalization of the primordial 
potential power spectrum k^P^o, and the angular size of 
the sound horizon at the last scattering surface 9 S . 

There arc three major difference between our choice of 
parameters and that of SE03. First, we do not include 
the matter fraction S l m because at fixed 9 S and ui m it af- 
fects neither the CMB observables nor the matter power 
spectrum. Second, we use lnG& instead of InG because 
the galaxy clustering bias is degenerate with the growth 
function. And third, we include additional parameters 
In a z (stage I) and 5z (stage II) in each redshift bin. We 
also replace -Dcmb in SE03 with 8 S for convenience, which 
is not critical. 

As the cosmic density field evolves, nonlinear effects 
become important on larger and larger scales. To pre- 
vent these effects from contaminating the measurements 
of baryon oscillations, we only use the Fourier modes 
at k < fcmax in our analysis, where fc max ranges from 
0.12 to 0.53 /iMpc -1 between z = 0.31 and 2.66 (see 
Table 0. Evidence from the 2-degree Field Galaxy 
Redshift Survey and the Sloan Digital Sky Survey sug- 
gests that the galaxy clustering bias at z ~ 0.1, despite 
some dependence on luminosity, is scale-ind ependent on 
scales k less than a few tent hs /iMpc" 1 fe.g. lVerde et al.l 
120021 iTegmark et"aD 12004ft . Hence, the upper bound 
in wavenumbcr also ensures that we can model the 
galaxy bias as a time-evolving but scale-independent 
quantity. For simplicity, our fiducial model has a bias 
b(z) = 1 + 0.84 z for all galaxies. It is straightforward 
to model the bias for each sub-sample of galaxies sep- 
arately. In fact, it is beneficial to divide the galaxies 
into sub-samples, especially at low redshift where galaxy 
number density is high, because one can control pho- 
tometric redshift errors much better for a homoge neous 
sample of galaxies (e.e. IPadmanabhan et al.|l2005jl . 

We assume that the primordial potential power 
spectrum fc 3 P|,(fc) = fc^P^fc/fco)" 5-1 - Our 
fiducial model is a low density and flat cold 
dark matter universe with a cosmological con- 
stant. It has ( uj m , w b , Zr e i ,y p ,n s ,k^P^ ,9 s ) = 
(0.146, 0.021, 6.3, 0.24, 1.0, 6.4 x 10" 11 , 0.60 deg) and the 



reduced Hubble constant h = 0.655. 

We calculate the Fisher matrix using the approxima- 
tion l|TegmarUll997t) 



F] 



i ainP(k f )ainP(k f ) 

j 2 dpi dpj ° f 

where pi is the ith parameter, 



dk f 

(2*0* 



Kff(k) = J 



i(r)P g (k) 



dr 



(4) 



(5) 



n(r)P g (k) + 

is th e effective volume of the survey (iFeldman et alJ 
Il994|) . and n(r) is the galaxy number density. The 
marginalized error of the ith parameter is given by [(F + 

Pp)" 1 ]-/ 2 . Here, F P is for any prior information. If we 
already know, prior to examining the data, the value of 
parameter pi to within ±ap(j>i) then (Fp)a = l/<7p(pj). 

In the second stage, our "observables" are ui m , 9 S , D's, 
and iJ's with their covariance matrix as calculated from 
the stage I Fisher matrix, marginalizing over all the other 
parameters. From these observables, and their covari- 
ance matrix, we forecast how well one can determine the 
following parameters: the photometric redshift biases 5z, 
u> m , QdEj wo, and w a , where the dark energy EOS is 
parametrized as w(z) = Wo + w a [l — (1 + z) -1 ]. To ex- 
plicitly define the redshift bias parameters we write the 
probability that a galaxy with photometric redshift es- 
timate z p has a true redshift in between z and z + dz 



dP(z) 



2tt a. 



exp 



(z + Sz - z p f 



2er?. 



(6) 



2.3. Tests 



We test our procedures using the baseline surveys of 
SE03. In order to compare with their results, we do 
not include photometric redshift parameters and do not 

differentiate the photometric redshift term e z u' 
with respect to H in the Fisher matrix. We also adopt 
SE03 parametrization of the dark energy EOS in this 
test. The resulting uncertainties on D and H match 
those in SE03 within a factor of 1.2 (1.45) for the spec- 
troscopic (photometric) baseline survey. Our error fore- 
casts on dark energy parameters are 30%-50% larger 
than those in SE03 because of the larger errors we obtain 
for D and H. 

If we include the derivative of the photometric redshift 

_ 2 2-.2 / ttI 

term e " ll' with respect to H in the Fisher ma- 
trix, then the constraints on H and dark energy EOS 
parameters are nearly as good as those from the spectro- 
scopic baseline survey. This is because the shape of the 
observed galaxy power spectrum is exponentially sensi- 
tive to the Hubble parameter. Such sensitivity imposes a 
strong constraint on H and subsequently on dark energy 
EOS parameters. However, the strong constraint is actu- 
ally on the comoving length scale of the rms photometric 
redshift error ca z /H. In reality, the uncertainty of a z 
will degrade the constraint on H. In fact, if we also in- 
clude hia z in the parameter set, the constraints become 
as weak as those from the photometric baseline survey 
in the previous paragraph. Thus, we conclude that the 
photometric results in SE03 are equivalent to having an 
imprecise prior on a z ; i.e., letting a z float freely. 
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TABLE 1 

Marginalized Errors of Selected Stage I Parameters 
with Complete Self-Calibration from LSST 



z 


femax 
(h Mpc" 1 ) 


cr(lllD) a 


o-(lnH) a 


o(\nH) h 


<r(ln a z /H)^ c 


o-(ln/3) a 


a(ln Gfc) a 








(X100) 






0.31 


0.12 


2.9 


27 


1.1 


0.37 


24 


13 


0.55 


0.14 


1.7 


19 


1.0 


0.21 


17 


9.5 


0.84 


0.17 


0.75 


11 


1.0 


0.14 


13 


7.1 


1.18 


0.21 


0.57 


12 


1.0 


0.11 


11 


6.2 


1.59 


0.28 


0.43 


10 


1.0 


0.10 


9,1 


5.2 


2.08 


0.38 


0.39 


11 


1.0 


0.11 


10 


5.5 


2.66 


0.53 


0.43 


12 


1.0 


0.15 


13 


6.4 



a No prior is taken for the rms photometric redshift error or the Hubble parameter. The 
constraints are poorer at low redshift due to small volumes and low fc max and then are poorer 
at high redshift due to low galaxy number densities. 

b We assume txp(ln<T z ) = 0.01 for this column. The prior sets the precision of H, reduces 
uncertainties of /3 and Gb significantly, and renders little improvement to other parameters. 

c For this column, the parameter lna z in the Fisher matrix is replaced by \na z /H. With 
this parametrization, the constraints on other parameters change slightly. 



3. RESULTS FOR REDSHIFT PARAMETERS 

We assume that through some external calibration pro- 
cess we have estimates of both the distribution of errors 
about the mean redshift as well as the mean redshift. 
A central point of our paper is that one can use the 
clustering properties of the galaxies themselves to test 
those externally-derived mean redshifts. To explain how 
it works, let us assume that the distribution of errors is 
Gaussian with variance a 2 . The redshift errors lead to 
a suppression of the fluctuation power of radial modes 
by an amount exp(— c 2 a z k 2 /H 2 ). With a z known, one 
can thus determine H. Analyzing the power spectra of a 
sequence of galaxy groupings, each with mean redshifts, 
z m , we thus know H(z m — Sz) for each group of galaxies. 
From the standard ruler provided by the baryon oscil- 
lations, we can also determine D(z m — 5z). With some 
assumptions about the smoothness of H , and using the 
fact that D = cj dz/H, these values of D and H will 
not be consistent with each other for all possible values 
of z m — Sz. Thus we can determine these mean redshifts 
and check to see if they are consistent with the mean 
redshifts inferred from the external calibration. 

Even if we do not use the external calibration to tell 
us a z , from the galaxy power spectra we can still deter- 
mine the source distance rms o~ r = ca z /H to much bet- 
ter than 1% (see Table Fortunately, this is exactly 
the quantity we want for interpretation of weak lensing 
power spectra, since the weak lensing power spectrum for 
sources with mean redshift z depends on how the sources 
are distributed in real space, not redshift space. 

In this section we investigate quantitatively how well 
one can measure tr z , given H (and vice- versa), perform 
the redshift bias consistency test and determine a z /H. 
For specificity, we adopt for this section and the next 
a model of the proposed survey by LSST as our fidu- 
cial survey. We assume a survey area of 23, 000 square 
degrees with galaxy distribution 

3 -z/0.35 „„„„,;„-2 



n(z) = QAOz 2 



aremm 



(7) 



which corresponds to a projected number density of 55 
arcmin -2 . We use a z o = 0.04 as our fiducial rms photo- 
metric redshift error. The survey is divided into 7 red- 
shift bins from z = 0.2 to 3 each with roughly the same 



comoving radial extent. The mean redshifts of the red- 
shift bins are listed in Table □ 

For the CMB we ass u me P lanck as treated in 
iKanlinghat. Knox, fc Sonel (|2003tL Planck observes at 
9 different frequency channels from 30 to 850 GHz. To 
crudely model the effects of foregrounds we assume the 
temperature maps at 100, 143 and 217 GHz can be per- 
fectly cleaned of foregrounds by use of the other channels 
and for polarization we keep just the 143 and 217 GHz 
channels. We assume this can be done over 80% of the 
sky up to / m ax = 2000 for temperature and l max = 2500 
for polarization. 

3.1. rms Photometric Redshift Error and H 

We first implement the full procedure described in Sec- 
tion §|2 with a prior on the Hubble parameter ap(hiH). 
For a simple demonstration, we apply the same prior 
to all redshift bins. The results for four of the red- 
shift bins are shown in Fig. ^ The other three are sup- 
pressed for clarity. As expected, the prior on H sets 
directly the uncertainty of the rms photometric redshift 
error, i.e. cr(lncr 2 ) = crp(lnff), over a wide range. One 
is limited by the uncertainty of the combination <j z /H 
when the prior crp(lnif) is stronger than the constraint 
<r[ln(cr 2 /iJ)], so that erp(lniJ) can no longer improve the 
precision of a z . On the other hand, cr(lner z ) does not de- 
grade further when the prior is worse than what can be 
determined from the survey itself. Hence, the a(hia z )- 
crp(lnff) curves in Fig. flatten on both sides. 

A prior on the Hubble parameter could potentially 
come from a spectroscopic galaxy survey such as that 
planned with the Kilo- Aperture Optical Spectrograph 4 . 
KAOS will be able to survey enough volume to determine 
the Hubble parameter at z ~ 1 and 3 to a few percent 
from baryon oscillations (SE03). 

On the other hand, we may take a prior on a z from 
a spectroscopic calibration of the photometric redshifts. 
We have thus also calculated the effect of a prior on the 
rms photometric redshift error ap(hia z ) on our ability 
to determine H. Again, the same prior was taken for the 
seven redshift bins. The results are almost identical to 

4 See |http://www.noao.edu7 kaos/ 
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Fig. 1. — Fractional errors on the rms photometric redshift er- 
ror cr(lncr 2 ) as a function of the prior on the Hubble parameter 
ap(lnff) in four of the seven redshift bins assuming the fiducial 
LSST survey. On the left side, the prior is better than the con- 
straint of the combination a z /H (see Table [Tt . so that it can no 
longer improve the precision of a z . Whereas, on the right side, 
the prior is not useful because the survey constrains H better than 
the prior, i.e. it reaches self-calibration. In the middle, one has 
er(ln (Tz) = crp(lnH) as expected. The uncertainty of H as a func- 
tion of the prior on a z look identical to this figure except the x-axis 
label becomes op (In a z ) and the y-axis label <r(ln //). 



those in Fig. Q] but with x-axis label becoming ap(lna z ) 
and the y-axis label a(\nH), because a z and H are de- 



generate in the photometric suppression e 



3.2. Photometric Redshift Bias 

In stage II of the Fisher matrix analysis, we include 
redshift biases Sz. These redshift biases will cause in- 
consistencies between the angular diameter distance and 
Hubble parameter in a given cosmological model, given 
some assumptions about the smoothness of H. In our 
analysis, those smoothness assumptions are provided 
implicitly by modeling the dark energy density as a 
smoothly varying function, controlled by EOS param- 
eters wo and iu a . 

The resulting constraints on the photometric redshift 
biases are shown in Fig. [21 as a function of the input 
prior crp(ln<7 z ), which is the same in all redshift bins. 
With crp(lncr z ) ~ 0.01 the biases are determined to bet- 
ter than 0.01. For t he 4,000 s q. degree cosmic shear 
survey considered by iMa et all 



(|2005D, these bias con- 
straints are sufficient to keep the degradation in wq and 
w a errors at the 20% and 50% levels respectively. How- 
ever, because the larger LSST cosmic shear survey can 
achieve stronger constraints on Wq and w a (given perfect 
redshift information), the degradation in wq and w a in 
this case is worse: factors of 2.5 and 3 respectively 5 . This 
degradation is much better though than the factor of 10 
degradation in wo and w a that occur without any prior 
on the redshift biases. 

At low <7p(lncr z ) the quality of the constraint as a func- 
tion of z peaks at intermediate rcdshifts, due to the simi- 
lar behavior for the constraints on D and H (see Tablc^). 

5 Zhaoming Ma, private communication. 




lO" 3 0.01 0.1 1 

(7 p (ln ct z ) 

Fig. 2. — Forecasted errors on Sz given varying priors on cr z . 
At low up^no^) the quality of the constraint as a function of z 
peaks at intermediate rcdshifts, due to the similar behavior for the 
constraints on D and H (see Table QJ. 



4. DARK ENERGY CONSTRAINTS FROM LSST BARYON 
ACOUSTIC OSCILLATIONS 

While the constraints on o~ z and 5z from galaxy power 
spectra discussed in Section [3] will be useful for extract- 
ing dark energy constraints from cosmic shear data, we 
also explore how well dark energy information can be 
extracted from the galaxy power spectra themselves. 

In Figures [31 and 0]we see contours of constant fore- 
casted errors in wq and w a as the redshift priors are 
varied. The same priors are applied to all redshift bins 
for simplicity. As expected, the galaxy power spectrum 
alone provides for only a very weak self-calibration; as 
the prior on rms redshift errors is relaxed to 100% and 
that on redshift biases to 0.1, the errors degrade by fac- 
tors of 20 to 40. Also as expected, we do see though that 
with a strong prior on cr z , one does not need a prior on 
the 5z parameters. 

From Fisher matrix calculations like those in Section[21 
we find that the Hubble parameter does not significantly 
improve the constraints on dark energy EOS parame- 
ters unless it is determined with precision comparable to 
or better than those of the angular diameter distance. 
This means that one must achieve cr(lniT) < 1%, which 
is possible, as discussed in Section 13.11 with a prior of 
op(lna z ) < 1%. 

We show LSST baryon oscillation constraints on wo 
and w a in Fig. 0for two cases: [<r z , crp (In <j z ), ap(Sz)} = 
(0.04,0.01,0.01) (dashed line) and (0.08,0.005,0.08) 
(solid line). By reducing crp (In a z ) from 0.01 to 0.005 one 
tightens the constraints considerably, even though the 
rms photometric redshift error is increased. This exam- 
ple demonstrates that a large rms photometric redshift 
error can be tolerated as long as it is known accurately. 

We also show in Fig. 03 the constraints from the power 
spectra of LSST tomographic cosmic shear maps (dot- 
ted line) and 2000 supernovae as might be observed 
by the Joint Dark Energy M ission (dash-dotted line, 
IKnox. Albrecht. fc Songl 12005(1 . Note that these weak 
lensing error forecasts are done in the limit of perfect 
prior knowledge of the redshift error distribution param- 
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Fig. 3. — Contours of constant forecasted errors in too estimated 
from galaxy clustering in the LSST survey as priors on Sz and a z 
are varied. The same priors are applied to all redshift bins for 
simplicity. 
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eters. If we were to use the priors on these parameters 
adopted for the baryon oscillation calculations, the weak 
lensing constraints would become slightly looser than the 
baryon oscillation constraints. 

One may attempt to combine all the three constraints 
together to achieve higher precision. However, further 
work is needed to carefully account for the correlation 
between baryon oscillations and weak lensing statistics, 
because they overlap in b oth observational data and un- 
derlying density field (e.g. lHu fc Jainll200-|) . 

As mentioned earlier, we cannot yet be completely sure 
that the sub-bin evolution of galaxy bias can be con- 
trolled well enough to make an accurate determination 
of W (co z k\\ / H) , the suppression of power in the radial 
direction. If we are not able to extract the information 
about the shape of the window function and, therefore, 
a z /H, then even with a prior on the redshift error dis- 
tribution from external calibration our determination of 
H(z) will be degraded as discussed in Section UTTl We 
show in Fig. how the error contours degrade in the 
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Fig. 5. — Marginalized la error contours in the wg-w a plane. 
The contours for LSST baryon acoustic oscillations are shown 
for [<Tg,<TpQn<r z ) t crp(6z)] = (0.04,0.01,0.01) (dashed line) and 
(0.08,0.005,0.08) (solid line). For comparison, we include the er- 
ror ellipses from the power spectra of LSST tomographic cosmic 
shear maps (dotted line) and 2000 supernovac as might be observed 
by th e Joint Dark Energy Mission (dash-dotted line, IKnox et alj 
120051) . Note that the weak lensing constraints, unlike the baryon 
oscillation constraints, do not included the effects of redshift error 
distribution uncertainties. 




Fig. 6. — As in Fig. |3 but the two larger error contours are 
obtained without differentiating the photometric suppression term 
cxp(— c 2 a z k 2 /H 2 ) with respect to the Hubble parameter H in the 

Fisher matrix in stage I. The rms photometric redshift error is 
assumed to be a z = 0.04(1 + z), and it is not included as pa- 
rameters. In stage II, the same prior is applied to the photomet- 
ric redshift bias in each redshift bin. The solid contour corre- 
sponds to ap(8z) = 0.01, and the dotted contour op(<5z) = 0.001. 
The constraints do not improve for ap(Sz) < 0.001 or degrade 
much for crp(Sz) > 0.008. The dashed contour is from Fig. 151 for 
[a z , a P (In a z ),a P {Sz)} = (0.04,0.01,0.01). 



limit of no information extracted from the shape of the 
radial suppression. Specifically, for this calculation we 
do not differentiate the photometric suppression term 

e ' z ii ' with respect to H in stage I, and the parame- 
ters In a z are removed. In other words, this is the calcula- 
tion in the SE03 limit as discussed earlier. Unlike in SE03 
we include photometric redshift biases 5z in stage II to 
explore their effect on dark energy constraints. The error 
contour with [a z , a P (In a z ), a P (6z)} = (0.04,0.01,0.01) 
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(dashed line) from Fig. [3] is added for comparison. 

One sees in Fig. |H1 that the constraints from baryon 
oscillations in photometric redshift surveys are tight- 
ened tremendously by fully utilizing our knowledge of 
the photometric suppression to the galaxy power spec- 

_ 2 2 i 2 / tt2 

trum. Without extracting information from e = » ' , 
we obtain [cr(f2 DE ), <r{w ), a-(w a )] = (0.046, 0.33, 0.77) for 
a P {Sz) = 0.001 (dotted line) and (0.084,0.60,1.4) for 
ap(8z) = 0.01 (solid line). These constraints are as loose 
as those with very weak priors on In a z in Figs. |3] and 
0| They do not improve for ap(5z) < 0.001 or degrade 
much for a P (Sz) > 0.008. 

5. DISCUSSION AND CONCLUSIONS 

The statistical properties of weak lensing maps are sen- 
sitive to the distances to the sources, and not to the red- 
shifts of the sources. Therefore, even without redshift 
information, weak lensing observations can be used to 
determine the source distance distribution. The redshift 
information is needed only so the source distance distri- 
bution can inform us about the distancc-rcdshift relation. 
However, this distance-redshift relation is exactly what 
we want, because of its sensitivity to dark energy, and 
therefore the redshift information is actually terribly im- 
portant. 

We have pointed out that the three-dimensional power 
spectra of the source galaxies, because of their sensitiv- 
ity to redshift errors, can be used to inform us about 
the source redshift distribution. While the galaxy power 
spectra, on their own, cannot give us all we need to know 
for interpretation of tomographic cosmic shear data, they 
can be used to supplement the redshift information avail- 
able from external calibration of the photometric red- 
shifts. 

We have demonstrated the supplemental information 
available in two different ways. First, if one uses prior in- 
formation (perhaps from external calibration) about the 
variance of the redshift errors in each of the redshift bins, 
then the mean redshifts of each of those redshift bins can 
be recovered. They can be recovered with sufficient pre- 
cision to greatly limit the degradation in the errors on 
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dark energy EOS parameters reconstructed from cosmic 
shear data. Second, without any prior on the redshift er- 
ror variance, the three-dimensional galaxy power spectra 
can be used to determine very precisely the source dis- 
tance distribution which with the assumption of a nor- 
mal distribution is specified by the combination a z /H. 
The source distance distribution will be very helpful for 
interpretation of cosmic shear data. Although, for de- 
termining the distance-redshift relation (and thereby the 
dark energy EOS parameters) we will still require mean 
redshifts for each of the redshift bins, in this case from 
some external calibration. 

The galaxy power spectra, combined with external cal- 
ibration of the source redshift distributions, are also po- 
tentially powerful probes of dark energy. If the rms pho- 
tometric redshift error is determined to better than 1% 
through some external calibration of an unbiased sub- 
sample, and one can control sub redshift bin bias evolu- 
tion then one can significantly reduce errors on the Hub- 
ble parameter and, consequently, put tight constraints on 
the dark energy EOS parameters using baryon acoustic 
oscillations. For instance, with priors erp(lner 2 ) = 0.005 
and ap(Sz) < 0.004, one can achieve cj(£!de) = 0.015, 
o-(wo) = 0.095, and cr(w a ) = 0.24 with LSST. 

Achieving, in practice, the constraints on dark energy 
from baryon oscillations and cosmic shear, that are pos- 
sible in principle, will require very accurate calibrations 
of photometric redshifts for very faint galaxies. This is 
a significant challenge for the observational community. 
The consistency tests we have described here may play 
an important role in meeting that challenge. 
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